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ABSTRACT 

In this letter we present a study of the size luminosity relation of 475 early-type galaxies in the 
Virgo Cluster with Sloan Digital Sky Survey imaging data. The analysis of our homogeneous, model- 
independent data set reveals that giant and dwarf early-type galaxies do not form one common se- 
quence in this relation. The dwarfs seem to show weak or no dependence on luminosity, and do not 
fall on the extension of the rather steep relation of the giants. Under the assumption that the light 
profile shape varies continuously with magnitude, a curved relation of size and magnitude would be 
expected. While the galaxies do roughly follow this trend overall, we find that the dwarf galaxies are 
significantly larger and the low-luminosity giants are significantly smaller than what is predicted. We 
come to the conclusion that in this scaling relation there is not one common sequence from dwarfs to 
giants, but a dichotomy which can not be resolved by varying profile shapes. The comparison of our 
data to a semi-analytic model supports the idea of a physical origin of this dichotomy. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: dwarf — galaxies: fundamental 
parameters — galaxies: clusters: individual: (Virgo Cluster) 



1. INTRODUCTION 

Scaling relations have ever been an important tool not 
only to study galaxy properties but also to link those 
properties to their formation and evolution. Since early- 
type galaxies are the most numerous galaxy type in clus- 
ter environments, they play an outstanding role in un- 
derstanding galaxy clusters and therefore structure for- 
mation in general. But still today it remains an open 
question whether the giant early-type galaxies and their 
lighter counterparts share the same origin and forma- 
tion mechanisms. Of the morphological scaling rela- 
tions, so far mostly the relation betw een surface bright - 
ness and size ("Kormendy relation", lKormendvlll985D . 
' and between surface brig htness and luminosity (e.g. 
\ iBinggeli fc CameronI Il991f ) were studied to tackle this 
question. Combin ed with velocity dispersi on, the Faber- 
Jackson relation (jFaber fc JacksonI I1976D and the ex- 
\ tension to the Fundamen tal Plane (jPressler et al.|[T987l : 
. iDiorgo vski fc Davislll987l ) were studied, mostly for giant 
early types. But even with the now available facilities 
velocity dispersions are still rare for dwarf galaxies. 

Until the early 1990's sizes of near by early-typ e galax - 
ies were stud i ed fo r example by iKormendvl (|1977D . 
iGuzman e t al. (1993), and for the Virgo Cluster in par- 
ticular bym inggcli fc Cameron (1991) for dwarfs and by 
ICaon et all (jl993l ) for giants. An often cited source of a 
homogenous data set of sizes for dy namically hot s ystern s 
over the whole luminosity range is iBender et afl i|l992f ). 
The conclusion of that time was that giant and dwarf 
early-type galaxies show a distinct size distribution, and 
that the dwarfs show less change of size with luminosity 
than the giants. This together with other scaling rela- 
tions was interpreted as evidence for a different origin of 
dwarf and giant early- type galaxies. 

Towards the turn of the millenium, however, it be- 
came more widely realized that the light profile shapes 
of early types vary continuously with luminosity. Nei- 
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ther do dwarf galaxies simply follow exponential pro- 
files, nor do all giants exhibit de Vaucouleurs pro- 
files. Instead, all early ty pes are well described by 
the generahzed S ersic profile (ISersiclll963l) with different 
Sersi c indices n ([Young fc Currid 119941 : iFerrarese et al.l 
l2006f ). Several authors reasoned that the scaling re- 
lations naturally follow what is predicted by n chang- 
ing linearly with magnitude, and th at all these galax- 
ies can indeed be of the sar ne kind (iJerien fc Binggelil 
19971: iB inggch fc Jerienlll99l iGraham fc Guzm an 200f 
Gavazzi et al.l l2005!). While this latter in terpretation is, 
for example, not shared by Bos cUi et al.l (|2QQS,), who ar- 
gue that dwarf and giant early-type galaxies have dif- 
ferent origins, these authors do agree that they can be 
seen as one structural family with a gradual variation 
of n with luminosity. Recently, iKormen dv et al.l ()2008l ) 
conclude that dwarfs and giants are structurally distinct, 
notwithstanding the Sersic continuum, which they con- 
sider insensitive to the physics dividing the two. In our 
analysis below we quantitatively pin down this structural 
distinction by taking the variation of profile shapes ex- 
plicitly into account. 

2. SAMPLE SELECTION AND IMAGING DATA 

Our sample is selected based on the Virgo Cluster 
Catalog (VCC: iBinggeh et al.|[l985l) . Only certain clus- 
ter members with rriB < 18.0 mag are taken into ac- 
count, which is the same magnitude limit up to which 
the VCC was found to be complete. This translates into 
Mb < —13.0 mag with our ado pted distance inodulu s of 
m-M=31.0 mag (d=15.85 Mpc. lGraham et al.lfT999( ). 

Galaxies listed as "SO:", "E/SO", "SO/Sa", and 
"SBO/SBa" are taken as SO, and one SO (VCC1902) is 
excluded, since it shows clear spiral arm structure. For 
the dwarfs, we selected galaxies classified as dE, dSO, 
and "dE:", whereas "dE/Im" as well as possible irregu- 
lars b ased on visual inspection are excluded (Lisker et al] 
l2007t l. We excluded 13 dwarfs where the Petrosian aper- 
ture (see below) could not be obtained, as well as 13 
dwarfs and four giants that are too strongly contam- 
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inated by the light of close neighbour objects. This 
leads to a working sample of 475 galaxies: 397 early- 
type dwarfs ("d Es"), 9 M32-ty pe candidates (as hsted 
in Table XIII of iBinggeh eFaLl ligSSl. and 67 E and SO 
galaxies. 

Th e Sloan Digital Sky Survey (SP SS) Data Release 
Five (|Adelman-McCarthv et al.|[2007[) covers all but six 
early- type dwarf galaxies of the VCC. The pixel scale 
of 0^596 corresponds to a physical size of 30 pc. For 
the analysis below, we use the r-band, which has the 
highest S/N. Since the quality of sky level subtrac- 
tion of the SDSS pipeline is insu fficient, we use sky - 
subtracted images as provided by iLisker et al.l (|2007f ). 
based on a careful subtraction method. The images 
were fiux-cali brated and corrected for Galactic extinc- 
tion (fSchlegel et al.lll998D . 

The image analysis is done largely in the same way as 
in iLisker et al.l (|2008D : For ea ch galaxy we d etermined 
a "Petrosian semimajor axis" (|Petrosianlll976l Op ), i.e., 
we use ellipses instead of ci rcles in the calcu lation of the 
Petrosian radius (see, e.g.. iLotz et al.ll2004D . The total 
fl-Ux in the r band was measured within a = 2ap , yield- 
ing a value for the half-light semimajor axis, ahi,r,uncorr- 
This Petrosian aperture still misses some fl ux, which is of 
particular relevance for the giant galaxies (Truiillo et al.l 
[2001). As improvement of Liskcr ct al. (2008), luminosi- 
ties and half li ght radii are correct ed for this missing flux 
according to Graham et al.l (|2005D . Axial ratio and po- 
sition angle were then determined through an isophotal 
flt at a = 2ahi.r- The effective radius is then given by 
r^jj — j,y/b/a with the axis ratio b/a. Additionally 

we fltted Sersic profiles. We omitted intensities at radii 
r < 2" in order to avoid seeing effects. 

Our data set is a very homogeneous set of parame- 
ters for galaxies in one cluster, based on data taken with 
the same instrument and reduced and analyzed with the 
same procedure. We point out that our derived radii are 
model independent.^ 

3. SIZES OF EARLY-TYPE GALAXIES 

In Fig. 1 we present the size luminosity diagram for our 
sample. At first glance the sequence from dwarf to giant 
early-type galaxies does not look very continuous: the 
giants follow a steep relation with a well-defined edge on 
the bright end of their distribution. The bunch of dwarfs 
apparently lie with a larger scatter around an effective 
radius oi r ^jj — 1 kpc, their sizes showing weak to no 

dependence on luminosity.^ 

Without applying the correction for the missing flux 
within the Petrosian aperture, the separation between 
the two sequences even widens (not shown) since the less 
compact objects of the same luminosity show less con- 
centration and therefore a smaller correction. 

It is important to notice that the different behaviour 
of dwarfs and giants does not depend on whether objects 

^ The corrections for the missing flux are based on the concentra- 
tion parameter of the galaxies; therefore they implicitely depend on 
the a ssumption of a (Sersic) hght profile model (cf. Graham et al.l 
[200l) . 

^ The almost constancy of the dwarfs is just a coincidence when 
taking the half light radius containing 50 % of the light. For ex- 
ample, for the radius containing 90 % of the light, the relation for 
dwarfs steepens. This can be understood when taking the gradual 
variation of profile shapes with luminosity (see below) into account. 
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Fig. 1. — Absolute magnitude in r versus logarithm of half light 
radius. dEs are shown with open triangles and E and SO with filled 
squares. M32-type candidates are drawn with open squares. 

with disk components are omitted. This can be seen 
from Fig. 2, where gray symbols indicate objects with 
disk components or disk-like structure. 

Instead of adopting one distance for all galaxies, the 
Virgo cluster can alternatively be described as a sub- 
structured system, with the different componen ts partly 
having different distances (jGavazzi et al.l [T999| ). Since 
this can affect the apparent size of a galaxy, we checked 
whether the distribution of galaxies within the size lu- 
minosity diagram correlates with projected position in 
the cluster. We do not find any such correlation. There- 
fore, while distance variations could explain the larger 
scatter of the sizes of early types in the Virgo Cluster as 
compared to other clusters (jBosel li ct al. 2008) , it cannot 
explain the observed dichotomy: [Boselli et aLi 's Fig. 6 in- 
dicates the apparent separation of dwarfs and giants and 
they use different distance moduli for different subparts 
of the Virgo cluster. 

The same impression c an also be obtained from 
other pr evious studies , e.g. 'Bing geli &: CameronI (|1991|). 
Fig, l b; iBender et all ([1992), fable 1: iKormendv et all 
(|2008f ). Fig. 37. It is, however, not as clearly seen in 
the compilation of sizes of elli ptical galaxies from s evera l 
different studies presented by iGraham fc WorlevI ()2008D 
(Fig. 10). In this more heterogeneous dataset, the rela- 
tive number of small low- luminosity giants as well as that 
of large bright dwarfs appears to be somewhat smaller. 

3.1. Varying profile shapes ? 

iGraham Sz GuzmanI ()2003f ) suggested that the appar- 
ent dichotomy between dwarfs and giants in scaling re- 
lations can be explained just by the fact that the pro- 
file shape of a galaxy scales with magnitude. They de- 
scribe the light profiles with Sersic profiles and show the 
effect of a linear relation between magnitude and log- 
arithm of the Sersic index n on the other scaling rela- 
tions. As a result, the dependence of effective radius 
on magnitude becomes stronger at higher luminosities 
a nd the brightest galaxies are naturally larger (Fig. 11 
in IGraham fc Worlevl[200l . 

For investigating whether our galaxies display the pre- 
dicted behavior, we fitted their azimuthally averaged 
light profiles with Sersic models (see Sect. The de- 
rived Sersic indices n and central surface brightnesses fj,o 
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Fig. 2. — Left panel: Same as F ig. 1, wit h prediction for Sersic profiles with the Sersic index n given as a function of luminosity. Now 
the sub-classifications according to ILisker~et al. (2007) are taken into account. The symbols mean: filled squares - E, gray stars - SO, open 
squares - M32 candidates, open triangles - dE,N, open pentagons - dE,nN and gray open circles for dwarf galaxies with probable disk-like 
structure (dE{di) or dE(bc)). Right panel: The residuals in log(r) about the predicted relation. The vertical histograms show the galaxy 
distributions in magnitude bins of 1.2 mag. Their zero levels are indicated by the thin gray lines. Thick gray lines measure the scatter 
around the mean (black marks), with the standard deviation indicated in black at the ends of the lines. 

the largest scatter. Moreover, the distribution around 
the mean changes from the Gaussian-like shape seen at 
the bright and faint end to a much broader, even double- 
peaked shape. For those two bins, a K-S test yields a 
probability for the null hypothesis that the two samples 
follow the same distribution of 2.11%, which means that 
the break is statistically significant. 

Only the bright end of the predicted curve shows strong 
sensitivity to the fitted relations. At intermediate lumi- 
nosities, it is stable against small changes of the fits. 
Moreover, a change of the curve would not reduce the 
significant difference between the bright dwarfs and low- 
luminosity giants. 

Thus our analysis shows two things. First, the residu- 
als do not resemble a quite large random scatter around 
the relation. And therefore the size luminosity relation 
can not be fully explained by varying profile shapes. Sec- 
ond, the abrupt change in the behavior of faint and bright 
galaxies is even emphasized through the above examina- 
tion, and this break is a real discontinuity of the sequence 
from lowest to highest luminosities. 

4. COMPARISON TO SEMI-ANALYTIC MODEL 



were then used to obtain linear fits to the iiq/Mj. and 
n/Mr relations, using a least squares fitting algorithm. 
For those fits we exclude systems with a (probable) disk 
component, namely galaxies classified as S O, dEs with 
disk fe atures, and dEs with blue centers (iLisker et al.l 
I2006bllah . This is to ensure that the light profiles can be 
well parametrized by Sersic profiles. Our fits together 
with equation (16) of iGraham fc WorlevI (|2008[ ) predict 
a non-linear sequence in the r^jj/Mr diagram. 

The predicted relation is shown together with the ob- 
served galaxies in the left panel of Fig. 2. With the visual 
guidance of the line, it appears more likely that the data 
points follow one common continuous relation. And the 
gross trend in the diagram can indeed be explained by 
varying profile shapes. However, at luminosities bright- 
wards of the transition between dwarfs and giants, a sub- 
stantial amount of galaxies fall below the relation, while 
faintwards most of the dwarfs lie above it. 

We note that M32-type candidates play a minor role 
just by their small number, and furthermore, they are 
fainter than the "compact giants" in question. Whether 
both of these are special enough to justify their own clas- 
sification is out of the scope of this letter, and will be 
investigated in a future study. 

To quantify the departure of the observed galaxies from 
the predicted relation we show in the right panel of Fig. 2 
the size residuals about the curve shown in the left panel. 
While at both ends of the luminosity range the observed 
galaxies fall onto the relation, it can now be seen even 
more clearly that towards intermediate luminosities, they 
depart more and more from the relation in opposite direc- 
tions. Brightwards of the transition region, most galaxies 
are only half the predicted size, while faintwards, many 
galaxies are larger than predicted by >50%. 

Furthermore, we divide the data into magnitude bins 
and investigate the distribution of galaxies in those bins 
with histograms (Fig. 2, right panel). The scatter and the 
shapes of the histograms confirm the coexistence of two 
separate relations. In the two bins in which the transition 
occurs, the scatter increases, and is even larger than for 
the faintest galaxies, for which one would naively expect 



The Numerical Galaxy Gatalog of iNagashima et al] 
(2005) is based on a hig h resoluti on iV-body simulation in 
a ACDM universe (iYahagill200,'Th . The dark-halo merger 
trees of the iV-body simulation are taken as input for a 
semi-analytic model of the physical processes governing 
galaxy formation and evolution (here a mo dified version 
of the Mitaka Model, INagashima fc Yoshilll2004f l. 

In particular, this model takes into account the dynam- 
ical response to starburst-induced gas removal after gas- 
rich mergers (also for cases intermediate between a purely 
baryonic cloud and a baryonic cloud fully supported b y 
surrounding dark matter as in Yoshii fc Arimotolll987l) . 
This process plays a crucial role for the sizes of early-type 
dwarf galaxies. Their gravitational well is shallower and 
thus they suffer a more substantial gas loss than giants. 
The subsequent variation of the potential results in an 
increase in size. If it is not taken into account the dwarf 
galaxies are modelled to be systematically smaller. 

We identify model galaxies as early-type galaxies if 
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Fig. 3. — Absolute B magnitude versus logarithm of half 
light radius for o b served (black) and model galaxies (gray) from 
INagashima et all II2005I ). Model galaxies with fainter surface 
brightnesses than a limit of (/^s) < 25.5 mag/arcsec^ (dashed line) 
are excluded. Contours were calculated using model galaxy abun- 
dances in bins of 0.75 mag and 0.2 dcx. The contour levels, relative 
to the lowest, are 2.5, 5, 8, 11, 15, 20, 35, 50 

they are bulge dominated (bulge to total ratio > 0.6). To 
compare our data with the model we transformed SPS S 
g magnitudes into B according to ISmith et al.l ()2002D . 
using the galaxies' g — r color measured within ahi,r- In 
Fig. 3 one can see that the model galaxies show a bi- 
modality similar to what we observe, with low galaxy 
density between the two regions. In those dwarfs which 
form by gas-rich major mergers a starburst follows and 
the dwarfs are enlarged by the dynamic response to the 
subsequent gas loss. This mechanism is not at work 
in gas-deficient mergers, and the resulting galaxies are 
smaller. 

Note that INagashima et al.l assume de Vaucouleurs 
profiles to calculate projected half-light radii from half- 
mass radii. For exponential profiles, which would 
be more appropriate for dwarfs, the model galax- 
ies would shift upwards i n the diagram by 0.11 dex 
(INagashima fc Yoshiil [20031 ). 

It is inter esting, that t h ey co mpared their model to 
the data of ^ Bender et al.l (|1992D and concluded that a 
substantial fraction of model dwarfs is too large to match 
the observed galaxies. With our dataset, a much better 
agreement is found (also see Ide Riicke et al|[2005l) . 
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5. SUMMARY AND DISCUSSION 

We studied the size luminosity relation of Virgo clus- 
ter early- type galaxies, based on model- independent size 
measurements from SDSS imaging data. We find that 
the dwarfs do not fall on the extension of the rather 
steep sequence of the giants. While the gross trend in the 
size luminosity relation can be explained by light profile 
shapes becoming steeper for more luminous galaxies, a 
closer look reveals that there is a clear discontinuity in 
t he behaviour of faint a nd brigh t galaxies. 

iDabringhausen et al.l ([2008') compare dynamical 
masses and projected half-light radii of elliptical 
galaxies, bulges of spiral galaxies, dwarf spheroidals, 
massive compact objects (MCOs), and globular clusters 
(GCs). These authors also find a bimodality in the 
size distribution of elliptical galaxies. The "compact 
low-mass ellipticals" and the giant ellipticals lie on 
one sequence down to the MCOs and GCs, while the 
ellipticals of intermediate brightness in the larger branch 
lie on one seque nce with the dwarf spheroi dals. The 
interpretation of lOkazaki fc Taniguchil (|2000D . based on 
calculations of galaxy interactions, is that those latter 
galaxies are mostly of tidal origin and that the "compact 
low-mass ellipticals" should be interpreted as the real 
counterparts of the more massive elliptical galaxies. 
This is an alternative interpretation to the one based on 
semi-analytic models (Sect. |4]), and represents further 
evidence for the physical relevance of the observed 
dichotomy. 

Our findings are a new piece to the puzzle of the con- 
nection between dwarf and giant early-type galaxies and 
may hint at different origins. A further analysis of the 
scaling relations and comparison with models seems a 
promising approach to bring more light to this question. 
In particular it will be interesting to investigate correla- 
tions with velocity dispersion, since it should be closely 
related to size evolution in the framework of dynamical 
response. We will pursue this idea in a forthcoming pa- 
per. 
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